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ABSTRACT 

We present a detailed multi-wavelength study (from rest-frame UV to far-lR) of narrow- 
band (NB) selected, star-forming (SF) Ha emitters (HAEs) at z ^ 2.23 taken from the High 
Redshift(Z) Emission Line Survey (HiZELS). We find that HAEs have similar SED-derived 
properties and colors to sBzK galaxies and probe a well-defined portion of the SF population 
at z ~ 2. This is not true for Lya emitters (LAEs), which are strongly biased towards blue, 
less massive galaxies (missing a significant percentage of the SE population). Combining our 
Ha observations with matched, existing Lya data we determine that the Lya escape fraction 
(/esc) is low (only ^ 4.5% of HAEs show Lya emission) and decreases with increasing dust 
attenuation, UV continuum slope, stellar mass, and star formation rate (SER). This suggests 
that Lya preferentially escapes from blue galaxies with low dust attenuation. However, a 
small population of red and massive LAEs is also present in agreement with previous works. 
This indicates that dust and Lya are not mutually exclusive. Using different and completely 
independent measures of the total SER we show that the Ha emission is an excellent tracer of 
star formation at z 2 with deviations typically lower than 0.3 dex for individual galaxies. We 
find that the slope and zero-point of the HAE main-sequence (MS) at z 2 strongly depend 
on the dust correction method used to recover SER, although they are consistent with previous 
works when similar assumptions are made. 
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1 INTRODUCTION 

A wide variety of multi-wavelength surveys indicate that the rate at 
which galaxies form stars has changed with cosmic time, increas¬ 
ing by about one order of magnitude from Z'-^O t o z ~ 1—2, when 
the cosmic star formation had a likely maximum llLillv et alj T996I: 


[Perez-Gorizalez^alj[^05l: lHoDkins & BeacomlUOPd : Karim etZl] 

201 ll : Sobral et al 2013h . A similar behavior is found for the spe¬ 


cific star-formation rate (sSFR, the ratio between s tar-formation 
rate (SFR) and stellar mass) at a given stellar mass llNo eske et afl 
I 2 OO 7 I: iDutton et alJl2010l: iMagdis et al.ll2010l: ISobral et al.ll2014h . 
However, there is still some debate about the behavior at z > 2 — 3, 
since there are several uncertain factors involved in the analysis, 
such as dust correction factors or the c ontribution of emission lines 
in the determination of stellar mass (iBouwens et al.l l2009l . 1201 2l : 
IStark et ^l2013l[Gonzalez et al.ll2oT^ 


* Herschel is an ESA space obseivatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic¬ 
ipation from NASA, 
t E-mail: ivanoteogomez@gmail.com 


In agreement with the evolution of the cosmic SFR, the molec¬ 
ular gas content of star-forming (SF) galaxies is also higher at z ~ 2 
than at lower redshifts, although the star-formatio n efficiency is 
claimed to not strongly depend on cosmic epoch jTacconi et alj 
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l20ld:[Daddi et alj201(]l:lMagdis et al.l2012all9) . The implication of 
this is that the increase of the cosmic SFR density with redshift is 
most likely driven by an increase in the molecular gas mass fraction 
of galaxies. A change in morphology is also found, from disk-like 
shapes in the local u niverse to clumpy, irregular, or compact mor¬ 
phologies at z 2 telrnegreen et alll2009t iMalhotra et al.l 1201 2l : 
ISwinbank et al.l 12012 J bl). Given these rapid changes in galaxy 
properties it is important to ensure uniform selection of galaxy sam¬ 
ples at different epochs. 

Several selection criteria have been traditionally applied to 
select SF galaxies at z ~ 2: (1) the Lyoc narrow-band (NB) tech¬ 
nique, that selects Lyoc emitters (LAEs) by sampling the redshifted 
Lyot emi ssion line with a c ombination of NB and broad-band fil¬ 
ters (e.g. lOuchi et'^l2008ll : (2) the Lyman break technique, that 
selects Lyman-break galaxies (LBGs) by using colors that sam¬ 
ple the redshifted Lyman break and inc ludes a rest-frame UV 
color to rule out l ow-redshift interlo pers dSteidel et alj|2003h : (3) 
the BzK criterion llDaddi et alj|2004ll . that selects SF BzK (sBzK) 
galaxies with a double color selection criterion involving opti¬ 
cal and near-IR (NIR) bands; (4) the BM/BX criteria, that select 
galaxies within 1 < z < 3 with different c ombinations of optical 
broad-band filters dAdelberger et alJl2004h . Other selection crite¬ 
ria are based on FIR /sub-mm or radio data dChapman et alJl^OOSl ; 
iRiechers et al.ll201^ , although these identify only the highest SFR 
systems. The NB technique has been also applied in th e NIR regime 


I 995 I; Moorwood et all 2000l; 

Kurk et al. 2004 

; Geach et al.l 

200 a 

Haves et al. 2010b 

; Lee et al. 

l2012l Anetal, 

12014; ISobral et alj 

2013l;lTadaki et al.l 

2013h. In this way it is possible to also use Ha 


to select and study SF galaxies all the way from the local Universe 
up to z ~ 2 (with Ha moving from the optical into the K band), 
with this bein g a much simpler, se lf-consistent and well-understood 
selection (see lSobral et al.ir2013l the only work so far where Ha 
emission has been traced from optical to K band in a single data¬ 
set/analysis). 

In order to understand the nature of SF galaxies at the peak 
of galaxy formation and the biases and incompletenesses of each 
selection criterion, a comparison of the physical properties of the 
galaxy samples selected by the different techniques is required. 
This study is important for all evolutionary conclusions based 
on given population of galaxies, e.g. mass-metallicit y relations, 
gas fr actions, morphologies, colors, dynamics, etc dStott et alj 
I2013al f9). Furthermore, it provides a way to interpret galaxy evo¬ 
lution studies based on Lya and Lyman-break techniques, the ones 
that can be applied at the highest redshifts. 

A comparison of the properti es of LAEs, LBGs, a nd sBzK 
galax i es has been already done jGrazian et'H! 2007 1 _ ^_et_^ 
20091; IPentericci et al.l 2010l: Lv et alJl201ll;lHaberzettl et alji2012l: 


Oteo et all 2014h . Haberzettl et al. ( 1201 2h found that the BzK cri¬ 
terion is useful to select galaxies at z ~ 2, but the samples are bi- 
ased towards massiv e SF galaxies and with red stellar populations. 
iGrazian et alj ( l2007h report that the sBzK criterion is efficient in 
finding SF galaxies at z 2 but is highly contaminated by passively 
evolving galaxies at red z — K^ colors. They also fou nd that the 
Lyman -break criterion misses dusty starburst systems. lOteo et akl 
( I 2 OI 4 I 1 found that most LBGs can be selected as sBzKs, but most 
of the sBzK do not meet the Lyman-break criterion since the lat¬ 
ter is biased towards blue and/or UV-bright galaxies. Furthermore, 
they found that sBzK galaxies are similar to SF galaxies solely se¬ 
lected by their photometric redshift and, therefore, represent an ad¬ 
equate population to study the bulk of SF galaxies at z ~ 2. How¬ 
ever, the sBzK criterion cannot be used to carry out evolutionary 


studies, unlike the Lyman-break, Lya, or Ha ones, due to the use 
of a given filter set for the galaxy selection. However, extensions 
to higher redshif t have been proposed with other broad-band filters 
dGuo et alj|2012ll . Additionally, sBzK galaxies do not represent by 
themselves a purely SFR-selected sample, but have a more compli¬ 
cated selection function. 

The main objective of this paper is studying the properties a 
sample of HAEs at z ~ 2.23 selected from th e High Redshift(Z) 
Emission Line Survey (HiZELS iGeach et alj 1200^ ; ISobral et alj 
l2009blll I 2 OI 2 II 2 OI 3 II 2 OI 4 I 1 . HiZELS uses a set of NB filters in 
NIR bands to look for emission-line galaxies up to z ~ 9. The study 
of HAEs will also allow to analyse the accuracy of the Ha emission 
as a proxy of SFR and the relation between SFR and stellar mass 
at z ~ 2. Combining the Ha observations with available, matched 
Lya data we will also study the Lya escape fraction and its relation 
with galaxy properties. 

This paper is organized as follows. In ^we present the data 
sets used, the selection of our sources, and how we analyse them. 
In l^we study the nature of HAEs and compare them with LAEs, 
LBGs, and sBzK galaxies to place HAEs into the context of SF 
galaxies at the peak of cosmic star formation. In 0we study the 
population of galaxies with both Lya and Ha emission and analyse 
the Lya escape fraction at z ~ 2. We examine in ^the accuracy 
of Ha emission as a tracer of star formation at z ~ 2 and in ® 
we explore the location of our galaxies in an SFR-mass plane and 
discuss in detail the uncertainties in the determination of the slope 
of the main-sequence (MS) of star formation at z ~ 2. Finally, the 
main conclusions of the work are summarized in 0 

Throughout this paper all stellar masses and SFRs reported 
are derived by assuming a Salpeter IMF. We assume a flat universe 
with (Q.„,,O.A, hri) = (0.3,0.7,0.7), and all magnitudes are listed in 
the AB system dOke & Gunnlll98^ . 


2 METHODOLOGY 
2.1 Data sets 


Due to the availability of multi-wavelength data and NB Lya and 
Ha observations over an overlap ping redshift range, we focus in 
this work on the COSMOS field Jscoville et alj|200% . In order to 
sample the SEDs and study the stellar populations of the galax¬ 
ies analysed i n this paper ( see g2.2l (. we take optical to NIR pho¬ 
tometry from lllbert et a D ( l20I^. mid-IR IRAC a nd MIPS data 
from the S-COSMOS survey 1 Sanders et al.ll200% . and Herschel 
PACS and SPIRE data from PEP and HerMES projects, respec¬ 
tively dLutz et al.l201 ll ; loiiver et alj2012h . High-qualit y photomet¬ 
ric red shifts for the studied galaxies are taken from |llbert_et_^ 


j[201^. GALEX dZamoiski et al.l2()0% and CHANDRA dElvis et ^ 
I 2 OO 9 I 1 data will be also used. 


2.2 Source selection 

The main objective of this work is the analysis of the properties of 
NB-selected HAEs at z ~ 2.23. In some Sections of this paper we 
will use a sample of LAEs, LBGs, and sBzK galaxies at z 2 to 
help understand the properties of HAEs place them into the context 
of the SF population at z ~ 2. In this Section we explain how all 

these galaxies were selected. _ 

_ The sa mples of HAEs and L AEs are taken from ISobral et alj 

d2013h and iNilsson et aid d2009tl . respectively. LAEs and HAEs 
were selected via the NB technique, with an NB filter centered at 
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3963 A (129 A width) for LAEs and at 2.121 ^ml (210 A width) for 
HAEs. In addition to the NB criterion, HAE selection requires iden¬ 
tification of the detected emission-line as Ha at z ~ 2.23 rather than 
other line emitters at different redshifts. This selection includes a 
BzK cut to remove low-redshift interlopers, a Lyman-hreak like cut 
to remove z ~ 3.3 [OIII] emitters, and inclusion of double and triple 
line emitters from the combination of all HiZELS NB filters. It 
should be pointed out that the sBzK criterion applied to HAEs does 
not produce the loss of galaxies with unusual colors, but it is used 
to increase the completeness of the sample. Also, it might be possi¬ 
ble that a small number of HAEs has been missed because of their 
extremely blue SEDs (similar to those for LBGs, see ID. However, 
this percentage is estimated to be very low due to the use of double 
and triple line detections, since blue HAEs would have very strong 
emission lines and low red dening, and thus b e detectable in [Oil] 
or [OIII] as well as Ha Isee lSobral et al.ll2013L for more details). 

Completeness analysis indicates that LAEs are 90% complete 
down to aLya flux of fi^ya ~ 6 x 10^*^ ergcm^^ s^* ( INilsson et alj 
l2009h and HAEs are 90% complete down to a n Ha flux of /na ~ 
5.6 X 10^’^ergcm^^s^* dSobral et alj 12013). We recall at this 
poi nt that the Lyg and Ha NB filters used in Nilsson et alJ ( l2009l) 
an d lSobraletaH l l2013h . respectively, select galaxies over an over¬ 
lapping redshift range. The Lya filter is broader than the Ha filter in 
the redshift space and, therefore, selects galaxies over a wider red- 
shift range, with the redshift range of HAEs being fully included 
within the redshift range of LAEs. The great advantage is that it 
will be possible to study galaxies with both Lya and Ha emission 
even if Lya has a velocity offset with respect Ha (see ©. 

It is im portant to point out t hat among the whole sample of 
18 7 LAEs of Nilsson et alj ( l2009h . only 118 have a counterpart in 
the lllbert et alj j2013t) catalog, the data set that we use for optical- 
to-NIR SED fits. This represents 63% of the sample. The non¬ 
detections are mainly a consequence of the blue nature of these 
LAEs, that are clearly detected in U, B, r, or i bands, but are very 
faint in z and redder bands. The non detection in the NIR is an in¬ 
dication of their low stellar mass, being less massive than HAEs 
and other SF galaxies detected in the NIR. The significant number 
of LAEs without NIR counterpart indicates that the Lya technique 
tends to select low-mass, blue galaxies. Accurate SED fits cannot 
be carried out for those faint LAEs due to the lack of near-IR in¬ 
formation that is essential for age, stellar mass, and redshift esti¬ 
mations. Stacking might be an alternative, but it has been reported 
that stacking in LAEs does not provid e reliable estimatio ns of the 
median properties of the population JVargas et alj|2014l) . There¬ 
fore, we have decided not to include these LAEs in the analysis. 
This might bias our results since we only include in the final sam- 
ple the most massiv e LAEs selected through the NB technique in 
INilsson et alj l l2009l) . We will indicate the implications of this in 
the relevant sections of the paper. 

Regarding H AEs, most of them have a counterpart in the 
lllbert et S] ( 1201 3l) catalog, with only 8% being undetected. The 
non detections are mainly due to the faintness of that small pop¬ 
ulation of HAEs in optical bands. At the same time this is an in¬ 
dication that the HAE selection is also able to identify very dusty 
sources that might be missed in optical-based studies. We do not 
include the previous 8% of faint HAEs in our analysis since their 
UV continuum cannot be well constrained. They represent a very 
low percentage of the total sample and, therefore, their exclusion 
is not expected to affect significantly the conclusions presented in 
this work. However, it should be noted that due to their red colors 
(e.g. high dust extinction) they might be a significant proportion of 
dust-extinguished HAEs. 


LBGs and BzK galaxies are taken from lOteo et ^ ( |2014|) . 
LBGs were selected with the classical dropout technique, where the 
NUV and U bands were used to sample the Lyman break at z ~ 2 
and a (/-V color was used to avoid contamination from low-redshift 
interlopers. Furthermore, non detection in the GALEX FUV band 
w as imposed. The B zK galaxies were selected using the criterion 
of iDaddi et aT] | |2004|) . that picks up both SF galaxies (sBzK galax¬ 
ies) and quiescent galaxies (pBzK galaxies). Since we are interested 
in galaxies dominated by star formation, we only consider sBzK 
galaxies for most of the analysis presented in this work, although 
pBzK galaxies will be used for a comparison in some discussions. 
While LAEs and HAEs have a narrow redshift distribution due to 
their selection with NB filters, LBGs and sBzK have redshifts span - 
ning typically 1.5 < z < 2.5 dPaddi et alj|2004l : lOteo et alj|2()14l) . 
Therefore, in order to carry out a fairer comparison with LAEs and 
HAEs at z 2.23, we additionally limit the photometric redshift of 
LBGs and sBzK galaxies to 2.0 < Zphot < 2.5. This redshift range 
has been selected to account for th e uncertainties of photometric 
redshift determinations at z 2.25 ( lllbert et al.ll2013h . To avoid a 
possible presence of any low-redshift interlopers in the HAE and 
LAE samples, we also limit their photometric redshifts to the same 
range: 2.0 < Zphot < 2.5. Again, this range is chosen to account for 
the Zphot uncertainties. 

We clean all samples from AGN contam ination by remov ing 
all sources with X-ray CHANDRA detections dElvis et '^120091) . It 
should be pointed out that the percentage of X-ray detections in 
our galaxies is low, less than 5% in all t he four samples studi ed. 
Additionally, we use GALEX photometry dZamoiski et alj|200% to 
clean our samples from low-redshift interlopers. At z ~ 2 the Ly¬ 
man break is redshifted out of the UV regime and, therefore, our 
galaxies cannot be detected in GALEX bands. 

After all these considerations, we end up with a sample of 373 
HAEs, 69 LAEs, 3751 LBGs, and 13194 sBzK galaxies. We note 
that out of all the samples, HAEs are the ones drawn from the small¬ 
est volume, followed by LAEs, LBGs, and sBzK galaxies. Thus, the 
number of galaxies in each sample is largely driven by the differ¬ 
ent volumes. The number density of HAEs is 4.8 x 10^'^Mpc^^ 
down to log(LHa) > 42.0 and the number density of LAEs is 
1.6 X lO^^Mpc^^ down to log(LLya) > 42.3. The number den- 
sity of LAEs obtain ed here is smaller than the value reported in 
INilsson et alJ d2009h . mostly due to our inclusion of the criterion 
to clean the sample from lower-redshift interlopers and because 
we only include in the sampl e galaxies detected in the NIR. Fur¬ 
thermore, the value reported in iNilsson et alJ ( l2009l) was calculated 
over 28% of the area covered by their observations, which in turn 
represents 2% of the entire area of the COSMOS field. This makes 
their calculation very uncertain due to the influence of cosmic vari¬ 
ance. The number density of LBGs is 3.1 X 10 “^Mpc and the 
number density of sBzK galaxies is the highest in our samples, 
1.1 X lO^^Mpc^^, due to the high number of sources selected. 
Also note that LBGs and sBzK have different, more uncertain SFR 
limits. 


2.3 Analysis 

In the rest-frame UV to NIR regime we analysed the nature 
of o ur selected galaxies via th e traditional SED-fitting technique 
with iBruzual & Charloj (l2003l) (herea fter B C03) templa tes. To this 
end we used LePhare dArnouts et al. I l 19991: lllbert et al.ll2006h . We 
included the emission lines in the stellar popula tion templates 
dSchaerer & de Barrel 20091 : Ide Barros et al.[|20l^ since we are 
working with SF galaxies and, in fact, LAEs and HAEs are selected 
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through their emission lines. The strength of the Lya emission is 
more uncertain than other emission lines due to its resonant nature. 
Therefore, we do not include U -band information in the fits. This 
does not significantly change the values of the SED-derived proper¬ 
ties, since the UV continuum is well sampled with the other filters. 
In this way, the filters used in t he fits are: Subaru Bj, Vj, i+, z+, 
and V ISTA Y, J, H, and K, Albert et alJbOQ^: iMcCracken eTaP 
l20I2h. IRAC data fo r the IRAC-detected galaxies are also included 
J Sanders et alj[20()7h . 

The BC03 templates used in this work were built by assuming 
an exponentially declining SFH and a fixed metallicity Z = 0.2Zq. 
We considered time scale tspjj values of 0.1, 0.2, 1.0, 2.0, 3.0, and 
5.0 Gyr. We chose a fixed value for metallicity because this pa- 
rameter tends to suff er from large uncertainties (see for example 
Ide Barros et alj|20l^ . For age we considered values ranging from 
10 Myr to 3.4 Gyr, the age of the Universe at the median redshift 
of our galaxies. Ages values were taken in steps of 10 Myr from 10 
to 100 Myr, in steps of 20 Myr from 100 to 200 Myr, in steps of 50 
Myr from 200 to 500 Myr, in steps of 100 Myr from 500 Myr to 1 
Gyr, and in steps of 0.2 Gyr from 1. 0 to 3.4 Gyr. Dust att enuation 
was included in the templates via the lCalzetti et alj ( l2000h law and 
parameterized through the color excess in the stellar continuum, 
Es{B — V), for which values ranging from 0 to 0.7 in steps of 0.05 
were considered. 

Once the templates are fitted, the rest-frame UV luminosity for 
each galaxy was obtained from its normal i zed b est-fitted template 
and converted to SFRuv via the lKennicud l ll998[) calibration. Note 
that SFRuv is not corrected for dust attenuation. The UV contin¬ 
uum slope was obtained by fitting a power-law function to the UV 
continuum of each best-fitted template. Figure [T] shows the SED 
fit results for six HAEs randomly selected from our sample. For a 
reference, the best-fitted templates with no inclusion of emission 
lines (fitted with LePhare as well) are also included. It can be seen 
that emission lines have a clear effect on the observed fluxes and 
the best-fitted rest-frame optical continuum emission. This is be¬ 
cause the strongest rest-frame optical emission lines are sampled 
with some broad-band filters used in the fits; Ha is within the 
band, [OIII]5007 within the H band, and [OII]3727 within the J 
band. While the rest-frame UV SEDs are similar in all the cases, 
fainter rest-frame optical continua are obtained when including the 
effect of emission lines. This translates into lower stellar masses 
and younger ages. 

Among the whole sample of SF HAEs, only nine are individu¬ 
ally detected in any of the Herschel bands. This represents a detec¬ 
tion rate of 3%. Interestingly, despite being very low, the percentage 
of Herschel detections is higher for HAEs than for LBGs (~ 0.7 %) 
or sBzK galaxies (~ 0.5 %). This indicates that the Ha selection 
can recover not only relatively dust-free, but also highly dust ob¬ 
scured sources. However, the comparison between the number of 
Herschel detections and physical properties is challenging since at 
z ~ 2 Herschel only selects the most extreme galaxies rather than 
normal SF galaxies. Source confusion is also a major problem when 
analysing the FIR emission of UV, optical, or NIR-selected galax¬ 
ies. We have attempted to identify source confusion by analysing 
the optical ACS images of the galaxies along with the MIPS-24 
fjm and VLA contours. As an example, see results in Figure for 
the three SPIRE-500 /rm-detected HAEs. The size of the images 
are 40” on each side, slightly larger than the SPIRE-500 fjm beam, 
the Herschel band with the largest PSF. In the three cases, the lo¬ 
cation of the MIPS detection is coincident with a radio emission, 
indicating that there is no significant contribution of nearby EIR- 
bright sources that might contaminate the fluxes in SPIRE bands. 


The contamination is even more unlikely in PACS bands, since their 
PSF is 2-A times smaller than SPIRE beams. As a sanity check for 
SPIRE-detected sources, we have re-done the FIR SED fits includ¬ 
ing only their less-likely contaminated PACS photometry. The val¬ 
ues obtained for the total IR luminosity, and hence for SFRir, are 
in agreement with those including SPIRE data within the uncertain¬ 
ties (~ O.I — 0.2 dex). It should be pointed out that there are two 
LAEs individually detected in Herschel but they are also detected 
in X-ray and, consequently, have a li kely AGN nature and are not 
considered in this work (see 32.2l and [Bongiovanni et al.ll2010h . 


The IR SEDs o f the EI R-detected galaxies a re fitted with 
ICharv & Elb^j200lh (CEOl) . lDale & Heloul l l2002l) . |Polletta et alj 
( l2007h . and Berta et al.l ( 1201 3h templates. As an example, we show 
in Figure [3] the IR SEDs of the three HAEs detected in SPIRE- 
500 /jm. The presence of the rest-frame 1.6 fjm stellar bump (sam¬ 
pled with the IRAC bands at the redshift of our galaxies) is clear 
in the three galaxies, indicating that their SEDs are dominated by 
star formation. All the different templates fit well the IR SEDs of 
the galaxies, with similar values of x? • We choose to report the 
results obtained with CEOl templates, as in many previous works 
in the literature. The best-fitted CEOl templates are integrated be¬ 
tween rest-frame 8 and 1000 ftm to derive total IR lu minosities, 
Ltr. T hese are then converted into SFR by using the iKennicutl 
( Il998h calibration. The total SFR is then obtained by assuming that 
all the light absorbed by dust in the UV is reemitted in the FIR: 
SFRtotai = SFRuv + SFRir. 


Since most galaxies are not individually detected in the FIR, 
we also performed stacking analysis in Herschel bands to study 
the FIR emission o f //eric/ieZ-undetected galaxies (see for exam¬ 
ple [ftareT^ |2oU for stacking analysis in HAEs at z ~ 1.47). 
One single band near the peak of the dust SED is enough to es¬ 
timate the total IR luminosity. We focused on the PACS-160 fim 
band for stacking due to its relatively small beam and employ the 
residual maps, as in many previous works . We stacked by using 
the publicly available IAS Stacking Library dBethermin et al.l2010l) 
and uncertainties in the stacked fluxes were derived by using boot¬ 
straps. For HAEs, we stacked in different bins of stellar mass (from 
log(M*/iH 0 ) = 9.5 to 11 in bins of 0.5 dex) and dust-corrected 
Ha-derived SER (SERua) (from log (SFRHa[M 0 yr^']) = 1.0 to 
2.0 in bins of 0.5). Those bins cover the whole range of values for 
those parameters. For sBzKs, where the number of sources is high 
and allows stacking over more stellar mass bins, we stacked from 
\og{Mif / Mq) = 9.8 to 11.6 in bins of 0.2 dex. 

We only detected stacked fluxes for HAEs in the bin corre¬ 
sponding to the highest stellar mass ( 10 . 5 ^log(M,/M 0 ) < 11, 
where 66 galaxies are included and the stacked flux is /igo/jm = 
1.4±0.3 mjy, corresponding to log(LiR/L 0 ) = 11.8±0.1) and 
highest SFRea (1-5 < log(SFRHa) < 2 . 0 , where 166 sources are 
included and the stacked flux is /leo/im = TO ±0.2 mJy, corre¬ 
sponding to log(LiR/L 0 ) = 11.3 ±0.1). This is consistent with 
more mass ive SF galaxi e s bein g more affected by d ust, in agree- 
ment with iGarn & Best! (l2010l) . ISobral et al.l ( l2012h or llbar et all 
1 I 2 OI 3 I) . although it could be also due to a pure scaling of the SED. 
We also stacked LAEs as a function of stellar mass (with the same 
bins as for HAEs), but no stacked detections are recovered probably 
due to the low number of sources producing poor statistic and also 
to the less dusty nature of LAEs (see also In sBzK galaxies we 
recover detections for 9.8 < log(M*/M 0 ) < 11.0 with more than 
1000 galaxies in each bin. The recovered stacked fluxes are summa¬ 
rized in Table l2.3l The stacked PACS-160 fim fluxes were converted 
into Lir by using single band extrapolations with CEOl templates 
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Figure 1. Examples of SED fit results for 6 randomly-selected HAEs from our sample. These fits are representative of the whole sample of HAEs studied in 
this work. The best-fitted BC03 templates including emission lines are shown in orange, and best-fitted templates with no emission lines are in green. Red open 
squares ai‘e the observed photometric data. Due to the uncertain contribution of the Lya emission, the U band information has not been included in the fits. 
Our HAEs (and also LAEs) are selected for having strong Ha (and Lya) emission. Therefore, it is expected that emission lines affect their observed fluxes. It 
can be seen in this figure that the best-fitted templates including emission lines have fainter rest-frame optical continua and, therefore, represent less massive 
galaxies. If templates with no emission lines had been used, the stellar masses would have been overestimated. The UV continuum (and hence SFRuv) of the 
best-fitted templates does not change significantly when emission lines are included. 
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Figure 2. Analysis of source confusion and contamination at FIR wavelengths from nearby sources in the three SPIRE-500 ;um-detected HAEs. The MIPS-24 
fjm (green), and VLA (yellow) contours are over-plotted in the ACS I-band images of the galaxies. HAEs are located in the center of each image. Images 
are 40” on each side, slightly larger than the SPIRE-500 ^rm beam, the Herschel band with the largest PSF. The VLA and MIPS detections, that have better 
spatial resolution than PACS and SPIRE, indicate that there is no significant contribution of possible nearby FIR-bright sources that might contaminate the 
HAE Herschel fluxes. The PSE of PACS bands is much smaller than the SPIRE PSEs, so contamination is even more unlikely. 


Table 1. Stacked PACS-160 /rm fluxes and their associated total IR lumi¬ 
nosities for our sample of sBzK galaxies at 2.0 < Zphot < 2.5. 


Stellar mass range 

Stacked /i 60 ;jm 

log (Lir/Z, 0 ) 

9.8 log(M*/M0) < 10.0 

0.53 + 0.05 mJy 

11.4 + 0.1 

10.0 C log (M./Mq) < 10.2 

0.80 + 0.06 mJy 

11.5 + 0.1 

10.2 ^ log (M./Mq) < 10.4 

1.00 + 0.09 mJy 

11.6 + 0.1 

10.4^ log (M,/M 0 ) < 10.6 

1.20 + 0.12 mJy 

11.7 + 0.1 

10.6 ^ log(M./M q) < 10.8 

1.23 + 0.12 mJy 

11.7 + 0.1 

10 . 8 ^ log (M,/M 0 ) < 11.0 

1.40 + 0.17 mJy 

11.8 + 0.1 


and errors were obtained from the flux uncertainties. SFRjr were 
obtained with the lKennicud l ll998h calibration. 


3 THE NATURE OE Ha EMITTERS AT Z ~ 2 

Figure|4]shows the SED-derived dust attenuation, stellar mass, UV 
continuum slope, and dust-corrected SFR of our HAEs at z ~ 2.23. 
It can be seen that HAEs can be either dust-poor (and have blue UV 
slopes) or dusty (and have red UV slopes) and have a wide range of 
stellar masses and dust-corrected SFR. This already indicates that 
our sample of SF HAEs, selected down to a fixed dust-uncorrected 
SFR, is not significantly biased towards any SED-derived prop¬ 
erty. The distributions of SED-derived properties of HAEs resem¬ 
ble those for sBzK galaxies, one of the classical populations tradi¬ 
tionally used to study galaxy properties at z ~ 2. However, HAEs 
seem to have slightly higher stellar masses and SFR on average, 
probably due to their selection based on SFR. The distributions for 
HAEs contrast to those found for LAEs, which have much lower 
dust attenuation and stellar mass and much bluer UV continuum 
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Figure 3. IR SEDs of the three SPIRE-500 /jm-detected HAEs as an illustration of the SED-fitting results for the Her schel-se\sc\sA g a laxies studied in this 
work. We plot the obs erved photo metry in K^ , the four IRAC bands, MIPS-24 pm, PACS, and SPIRE. The best fitted ICharv & Elb3 j200lh . |Polletta et*^ 
j2007t) . lDale & HeloJ <2003) . and I Berta et alj j2Q13h templates are also shown, with the color code indicated in the bottom-right legend. In the SED fits, a 
redshift of z = 2.23 has been assumed. The data provide a good sampling of the dust emission peak in HAEs at z ~ 2.23. Therefore, the integration of the 
best-fitted templates between 8 and 1000 fmv provides an accurate determination of their total IR luminositie s and, consequently, d ust-corrected SFR. No 
significant difference is found for Lir when using different templates. We choose to report results obtained with ICharv & Elb^ j20Qlh templates, as in many 
previous works. 


HAEs 
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Figure 4. From left to right: Dust attenuation, stellar mass, UV continuum slope, and dust-corrected total SFR of our selected HAEs (orange filled histograms). 
We include the distributions obtained for LAEs, LBGs, and sBzK galaxies. These properties were obtained by fitting BC03 templates to the observed multi¬ 
wavelength photometry (from B to IRAC bands, when available) of each galaxy. The BC03 templates were buill by assuming an exponentially declining 
SFH and a fixed metallicity of Z = O.2Z0. Emission lines were also included in the templates. Histograms have been normalized to their maxima in order 
to clarify the representations. The distributions indicate that HAEs and sBzK galaxies have similar stellar populations, both in the median values and in the 
range covered. HAEs have a very-well defined selection criterion and, therefore, represent an excellent sample for studies of star formation at z ~ 2. LAEs are 
significantly biased towards blue, less massive, and dust-poor galaxies. 


slopes on average. This indicates that Lya and Ha samples are 
formed, on average, by galaxies with different stellar populations 
(although the values obtained for LAEs are within the distributions 
found for HAEs). This suggests that the selection based on Lya is 
much more biased than the selection in Ha and that using Lya to 
select high-redshift galaxies might lead to the loss of a significant 
population at a given redshift (mostly the reddest and most massive 
galaxies). Figure |4] shows that the sample of LBGs is also slightly 
biased towards galaxies with lower dust attenuation and bluer UV 
continuum slope, although the effect is not as strong as it is for 
LAEs. 

Figure[2 shows the relation between color and stellar mass for 
our HAEs. We choose the Y — color because it matches at z ~ 2 
with the rest-frame u — r color traditionally used to define the blue 


cloud and the red se quence in the local Universe (see for example 
IStrateva et al.ll200ll) . Defining a clear difference between the blue 
cloud and red sequence at z 2 is challenging due to the low num¬ 
ber of passive galaxies populating the red sequence. However, we 
consider here that pBzK galaxies represent the prototype of passive 
galaxies populating the red sequence at z ~ 2. These galaxies are 
also represented in FigurejS] HAEs cover a wide range of colors and 
stellar masses and, thus, they represent a diverse population with a 
range of properties. Ha selection only misses a small population of 
the bluest and least massive galaxies at z ~ 2.25. Some SF HAEs 
have colors even similar to pBzK galaxies, as ha ppens for obscured 
galax ies with intense, obscured star formation dOteo et ^l2013lj . 
|2014|) . 

Interestingly, although most galaxies with Lya emission have 
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Figure 5. Relation between the Y — K color and stellar mass for our sample 
of HAEs at z ~ 2.25. We also represent the location of LAEs, LBGs, and 
sBzK galaxies selected in the same field. At the redshift of our galaxies, 
Y — K matches the rest-frame u — r color used to study the col or distribution 
and d efine the blue cloud and red sequence in local galaxies iStrateva et alj 
l200lh . For a reference we a lso plot the location of the pBzK quiescent pop¬ 
ulation toaddi et alj|2004l) and LAHAEs (see |4). We also represent the 
distributions ofY — K color and stellar masses for each kind of galaxy but 
not for LAHAEs due to their low number. It can be clearly seen that HAEs 
are well distributed across a wide range of colors and stellar masses as are 
sBzK galaxies. However, LAEs and LBGs are the galaxies with the bluest 
colors and lowest masses, indicating the biased nature of those selections. 
LBGs and sBzK galaxies are distributed over a larger area due to their wider 
photomeh'ic redshift distributions. Most LAEs are blue and less massive, but 
there is a small population of LAEs with red colors. Although the number 
of such red LAEs is not high, it indicates that Lya emission can also escape 
from dusty, red systems, as previously reported from individual detections 
in EIR wavelengths. 


low dust attenuation and blue colors, there is a population of 12 red 
LAEs with log [M^/Mq) > 10.25. Their stellar masses are as high 
as the most massive HAEs. All 12 red LAEs are detected in IRAC 
and there is no indication of power-law like MIR SED, so these red 
LAEs are SE galaxies rather than AGNs (note that AGN contami¬ 
nation in all our samples was avoided hy discarding galaxies with 
X-ray detection). This population represents a low fraction of the 
whole LAE sample, but it indicates that Lya emission can also es¬ 
cape from dusty, massive, and red galaxies, as previously shown 
via optica l colors (?) and suhmm or EIR detectio ns at different 
redshifts ( QiagmanetalJ^ 20051 : lOteo et all l20I2al f3: ICasev et"^ 
I2OI2I : ISandherg et al. 20I5I) . Lva surveys over larger areas would 
be needed to increase the samples of red and massive LAEs and 
study in detail how and why Lya can escape from dusty galaxies. 

It is also important to examine the overlap between HAEs and 
other populations of SE galaxies at z ~ 2.25. In other words: what 
is the fraction of galaxies of a given type which could have been 
selected/missed through any of the other criteria? This study has 
been already don e for LBGs, UV-selected, and sBzK galaxies in 
lOteo et ahl ( l2014h . In that work, it was concluded that most LBGs 
can he also selected through the sBzK criterion and only 25% of 
sBzK galaxies would have been selected as LBGs (mainly due to 


Figure 6. Location of our HAEs at z ~ 2.2 5 in the color-colo r diagram 
employed to look for BzK galaxies at z ~ 2 IPaddi et alJl2007h . We also 
represent the location of LAEs, LBGs and LAHAEs (see 0. The distribu¬ 
tions of B — z and z^K colors are also represented for each kind of galaxy 
but not for LAHAEs due to the low number of galaxies in that sample. It 
can be seen that most LAEs, LBGs, and HAEs can also be selected through 
the sBzK criterion. There is a small population of LBGs (and less numerous 
sample of LAEs) that are missed through the sBzK criterion, indicating that 
it misses the bluest and youngest galaxies at z ~ 2. 


the bias of the Lyman-break selection towards UV-bright galaxies). 
Consequently, sBzK galaxies are a better representation of the bulk 
of SE galaxies over 1.5 < z < 2.5 than LBGs. 

Figure represents the classical color-color diagram em¬ 
ployed to selec t high-redshift BzK galaxies, both SE and pas¬ 
sively evolving dPaddi et alJ 120041) . Studying the location of our 
selected HAEs in such a diagram gives information about the over¬ 
lap between those populations. Most HAEs can be also selected 
as sBzK galaxies. As commented in ^2.21 this was expected for 
HAEs since a BzK selection was applied in fSobral et alj 12013 to 
the NB-selected HAEs to avoid contamination from low-redshift 
interlopers. HAEs cover the same range of colors than sBzK galax¬ 
ies confirming that they are not strongly biased either towards 
dust obscured or dust-free objects. However, as suggested by Fig¬ 
ure |5] the Ha selection might miss the bluest galaxies (now in 
terms of their z — K color) at z 2.25, as it happens to sBzK 
galaxies. Actually, there is a sub-population of ’blue’ LBGs with 
BzK = (z — A") AB - (S — z)ab -0.2 and (B — z)ab < 1 that can¬ 
not be selected as sBzK galaxies. This is mainly due to their blue 
z — A" color for their B — z color. At 2 < z < 2.5, the z — K color 
samples the Balmer and 4000 A breaks and, therefore, is a proxy of 
the age of the galaxies: younger galaxies should have bluer z —A" 
colors. This is further confirmed by the shape of their optical SEDs 
and ages, since 80% of the blue LBGs are younger than 100 Myr 
(always under the assumption of exponentially declining SFH and 
O.2Z0 metallicity and with the caution of the degeneracies between 
age, SFH, dust attenuation, and metallicity) with nearly flat SEDs 
(when flux densities are expressed per frequency units). Finally, 
the colors of LAEs indicate that they are more biased towards bluer 
galaxies, in agreement with the results obtained from SED fits and 
Figure |5l 
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Figure 7. Median SEDs of the HAEs, LAEs, LBGs and sBzK galaxies stud¬ 
ied in this work. Bottom-right panel shows the median SEDs renormalized 
to unity at 2.5 fim. In order to explore the EW of the Lya emission in each 
sample we show the region around Lya in the upper-left panel. In this panel, 
the location of Lya emission is indicated with a vertical dashed line. This 
Figure shows that LBGs have a bright rest-frame UV continuum, whereas 
sBzK galaxies are the most attenuated on average. Except for the difference 
in the normalization factors, the median SED of HAEs and sBzK galaxies 
are similar, in agreement with the similarity in their stellar populations. The 
upper-left panel indicates that the Lya emission boosts the U -flux only in 
LAEs, indicating that LBGs, HAEs, and sBzK galaxies do not have strong 
Lya with high equivalent width, on average. 


A significant percentage of HAEs cannot be selected as LBGs 
due to their faint rest-frame UV continuum. This indicates that a 
significant population of galaxies with intense Ha emission, and 
thus actively forming stars, would be missed if only the drop-out 
criterion had been applied. This, along with the relatively high per¬ 
centage of sBzK galaxies that could be missed in the drop-out cri¬ 
terion used in this work, suggests that the Lyman-break technique 
employed in this work is biased towards UV-bright galaxies and, 
consequently, would not be adequate for studying the properties of 
a complete census of galaxies at the peak of star formation. Th is 
result is importan t for evolutionary studies jStott et al. I l2013afi 
ISobral et al.ll20l4) . 

To further analyse the properties of our HAEs we study now 
their median rest-frame UV-to-optical SED and compare it to the 
median SED of LAEs, LBGs, and sBzK galaxies. The median SEDs 
were huilt from the median fluxes in each photometric hand. Then, 
following the same procedure presented in 32.3l we fitted their ob¬ 
served fluxes with BC03 templates including emission lines. In 
these fits we did not include the U band photometry because it 
might be potentially contaminated by the uncertain strength of the 
Lyot emission, not included in the templates. The results are shown 
in Eigure |7] Apart from the normalization factor, the shape of the 
median SED of HAEs and sBtK galaxies is similar, supporting pre¬ 
vious results derived with SED-derived properties and observe col¬ 
ors. The SED of HAEs constrast with the median SED of LAEs, 
which is much bluer, similar to what it happens to LBGs. This re¬ 
inforced previous claims that Lya selection is strongly biased and 


might miss a significant percentage of the SF population at the peak 
of cosmic star formation. 

As commented above, the U band might be affected by the 
uncertain strength of the Lya emission. The top-left inset plot in 
Eigure Q shows the region of the SED around the Lya emission. 
It can be seen that the U -band magnitudes of HAEs, LBGs, and 
sBzK galaxies agree well with the best-fitted templates. Therefore, 
in these galaxies the equivalent width of the Lya emission is not 
high enough to alter significantly the U -band fluxes. This indicates 
that HAEs, LBGs, and sBzK galaxies do not have Lya emission 
with high equivalent width on average. This is not the case for 
LAEs. The median U -band flux of LAEs is brighter than the pre¬ 
dicted by the templates, suggesting intense Lya emission with high 
equivalent widths as expected by their selection. These results in¬ 
dicate that the Lya emission has high equivalent width only in a 
small sample of galaxies, reducing the chances of finding strong 
Lya emission in a general population of SF galaxies at z ~ 2. Our 
selected LBGs have higher rest-frame UV luminosities than the 
other galaxies analysed in this work. Their lack of strong Lya emis- 
sion, on average, is t hus in agreement with the results reported in 
ISchaerer et al.l ( 1201 ih . who found that, at a given redshift, the Lya 
emission is more common in galaxies with fainter UV magnitudes 
(see also lStark et alJl2oTT . 

It should be pointed out that the comparison between the dif¬ 
ferent samples presented in this work is very dependent on the 
depth of the surveys used to select them. However, COSMOS is one 
of the deepest sets of multi-wavelength data available. Therefore, it 
makes a good judgement of the biases that affect samples selected 
at z ~ 2 through the studied selection criteria for most deep extra- 
galactic surveys. One alternative to overcome this limitation would 
be usin g the lensing powe r of massive clusters of galaxies, as was 
done in lAlavi et al.l ( l2014h . allowing to detect galaxies much fainter 
than previous surveys at z ~ 2. However, the latter is not the tradi¬ 
tional way of selecting LAEs, sBzK, or HAEs as normally they are 
searched and analysed in well-known cosmological fields where 
a wealth of deep multi-wavelength observations are available for 
their analysis. 


4 MATCHED LYa AND Ha EMITTERS: THE LYa 
ESCAPE FRACTION AT Z ~ 2.23 

Of special interest is the analysis of the galaxies which exhibit both 
Lya and Ha emission (Lya-Ha emitters, LAHAEs), not only due 
to the low number of such sources reported so far at z ~ 2, but also 
beca use they allow to constrain the escap e fraction of Lya pho¬ 
tons dHaves et al.ll20103 : ISong et al.ll201^ . As indicated in 32.21 
we can study LAHA Es because the Lya and Ha NB filters used in 
iNilsson et alJ (l2009ll and lSobral et al.l (l2013h select galaxies within 
an overlapping redshift slice over the same region of the sky. The 
Lya observations select galaxies over the redshift range 2.206 < 
ZLya ^ 2.312. The Ha observations cover 2.216 < ZHa Z 2.248. 
The Lya observations cover a wider redshift range than the Ha 
observations. Therefore, if HAEs had strong Lya emission, they 
should have been detected in the Lya filter. 

In the overlapping region of the COSMOS field observed ^ 
the Lya and Ha NB filters there are 158 HAEs and 146 LAEo 

* We should point out that we include here the whole sample of LAEs, 
both detected and undetected in near-IR wavelengths. The reason is that 
we are interested here in the fraction of HAEs with Lya emission, not in 
the SED-derived properties of galaxies with Lya emission (for which only 
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However, there are only 7 galaxies in common between the two 
samples. Therefore, only 4.5% of HAEs have strong enough Lya 
emission to be selected as LAEs. This implies very low Ly g es¬ 
cape fraction at z ~ 2 and agrees with iHaves et aP OOlOa) . who 
found six LAHAEs in their sample, implying an average Lya es¬ 
cape fraction of ~ 5%. These results reinforce the different nature 
of the galaxies selected through the Ha and Lya NB techniques 
and highlight the low chance of finding galaxies with Lya emis¬ 
sion, at least at z ~ 2.23. It should be remarked that the Hayes et al. 
results are based on Lya observations about 10 times deeper than 
those u sed in this work. Similarly, the Ha survey in iHaves et alJ 
( l 2010 al) is about 2 times deeper than HiZELS. However, we cover 
in this work an area of the sky that is about 20 times higher than 
in IHaves et"^ ( l 2010 al) . Therefore, the results complement each 
other in different regimes of Lya and Ha brightness and area sur¬ 
veyed. We have explicitly indicated in Eigure the location of 
LAHAEs. Among the seven joint detections, six have blue B — z 
colors compatible with almost flat UV continuum, whereas one of 
them has a red SED. Apart from the red LAHAE with SED-derived 
£ 5(6 — L) ~ 0.5, the remaining LAHAEs have a median dust atten¬ 
uation of £j(B — L) ~ 0.15. Their ages range within 1 - 900 Myr. 

We derive the Lya escape fraction in our LAHAEs from the 
ratio between the intrinsic and observed Lya luminosities and as¬ 
suming case B recombination, so the intrinsic Lya emission can 
be obtained from the intrinsic (dust-corrected) Ha luminosity. In 
order to make a fair comparison with previous works, the Ha lumi¬ 
nosity is corrected for dust attenuation by using the SED-derived 
EfB-Vy. 


r _ _ ^obs {Eya.) _ 

~ 8.7 xLobs(^CX) X 100-4x£dB-V)xtr(XH„) 

In the previous equation we have assumed that Es(B — V) = Eg{B — 
V), as it has been reported in s everal previous works t o happen 
at hig h redshift (see for example TReddv & Steidel|[200l : lErb et alj 
l2006t[Roddv et alJl201Ci) . However, it should be remarked that this 
has not been proven accurately, and will represent one of the major 
sources of uncertainties in our derived Lya escape fraction. This 
uncertainty affects not only our results but also the results reported 
in most previous works at similar redshifts, since a relation between 
Es{B — V) and Eg{B — V) must be assumed. 

In Figure[ 8 ]we represent the Lya escape fraction as a function 
of the SED-derived dust attenuation, UV continuum slope, stel¬ 
lar mass, and dust-corrected SFRna for our LAHAEs at z ~ 2.23. 
For a comparison we also represent the nine LAHAEs reported in 
ISong et all i2014h with both Lya and Ha emission ( UV contin¬ 
uum slopes are not provided for individual galaxies in ISong et alj 
l2014h . To be consistent, we take the observed Lya and Ha fluxes 
and the SED-derived properties from ISong et ^ ( 12014h and then 
the Lya escape fraction is calculated with the previous equation. 
However, we note that the selection criteria are not the same in 
both samples. Here we use the classical NB technique to lo ok for 
LAEs and HAEs, whereas the LAEs in ISong et ^ | |2014|) were 
found through blind spectroscopy and then followed-up with near- 
IR spectroscop ic observations. A ctually, the Lya luminosities of 
the LAHAEs in ISong et al.l ( |2014|) are higher than for the LAHAEs 
presented in this work. Despite these differences, the Lya escape 
fractions derived in both works agree well. 

It can be seen from Eigure that the Lya escape fraction 


near-IR-detected LAEs should be included - see j|2.2F The seven LAHAEs 
have detection in near-IR wavelengths, so good SED fits can be carried out. 


decreases with increasing dus t attenuation, as previously reported 
at similar and lo wer redshifts dHaves et alj|20l4 lAtek et alj|20l4 
ISong et ^|2014|) . and also with increasing (redder) UV continuum 
slope (a proxy for dust attenuation) and higher stellar mass and 
SFR. We also include in Figure [ 8 ] the upper limits on the Lya es¬ 
cape fraction in those HAEs whose Lya is not detected. These up¬ 
per limits have been ca lculated assuming a limiting Lya EW of 80 
A dNilsson et alj|200^ . The upper limits indicate that the relations 
found for LAHAEs corresponds to the upper envelope of the actual 
correlations. However, upper limits also indicate that the highest 
Lya escape fraction seen at low dust attenuation and/or UV con¬ 
tinuum slopes are not seen at higher dust attenuation and/or UV 
continuum slopes. These results indicate that Lya emission prefer¬ 
entially escapes from blue, low-mass galaxies with low dust atten¬ 
uation (with the exception of a low percentage of massive, red and 
dusty LAEs as obtained in This is in agreement with the SED- 
derived properties and colors found for LAEs in previous sections. 
It should be noted that the number of LAHAEs at z ~ 2 reported so 
far is low and, consequently, matched Lya-Ha surveys over larger 
areas of the sky are needed to increase the number of LAHAEs and 
have more robust results with higher statistical significance. 


5 Ha SFR VS UV AND UV+IR: A CONSISTENT VIEW 


The Ha emission is an excellent tracer of instantaneous star for¬ 
mation and calibrations between the SFR and the Ha luminosity 
have been proposed in the literature (see for example iKennicud 
Il998l) . The Ha emission is affected by dust attenuation (although 
to a much lesser extent than, for example, the UV). Therefore, in 
order to derive the total SFR from Ha, accurate dust correction 
factors are needed. In the local universe, where a significant per¬ 
centage of galaxies can be detected in the FIR or their Ha and Hp 
emissions can be measured, the dust correction factors can be deter¬ 
mined with acceptable accuracy. However, this is much more chal¬ 
lenging for galaxies in the high-redshift Universe, complicating the 
determination of the total SFR. 

We examine in this section the robustness of the Ha emission 
as a tracer of SFR at z ~ 2.23 for our HAEs. We first compare the 
results obtained from Ha and UV estimates of the SFR when dust 
correction is not taken into consideration. This is shown on the left 
panel of Figure As it can be seen, the two tracers of star forma¬ 
tion give significantly different results, with the Ha-derived SFR 
being higher than the obtained from the rest-frame UV luminos¬ 
ity. The median values found are SFRna = 13.7±9.OM0yr^' and 
SFRuv = 7.4±6.6M0yr^^, where the uncertainties represent the 
interquartile ranges. 

In the right panel of Figure [9] we show the comparison be¬ 
tween Ha-derived and UV-derived SFR when dust correction is 
included. We correct the rest-frame UV luminosities by using the 
relation between the dust at tenuation and UV continuum slope de¬ 
rived in iHeinis et alj ( 1201 3l) . as in ^ The dust correction of the 
Ha luminosity has been derived by using the r elatio n between stel¬ 
lar mass and dust attenuation of I Gam & Best! ( l2010l) . that has been 
suggested to be valid for HAEs at least up to z ~ 1.5 dSo bral et af] 
I 201 I llbar etalJlioH: iDominviiez erai]l2013l: IPrice et alj 120141) . 


In the right panel of Figure |3 it can be seen that there is a very 
good agreement between the UV-derived and the Ha derived SFRs, 
despite these being calculated with completely different and in¬ 
dependent estimators at different wavelengths. The difference be¬ 
tween the two estimators is typically lower than 0.3 dex for in¬ 
dividual galaxies. Furthermore, both estimations agree quite well 
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Figure 8. Lya escape fraction as a function of dust attenuation {upper left), UV continuum slope {upper right), stellar mass {bottom left), and dust-comected 
SFRun, {bottom r ight) for our sample of LAHAEs (orange filled dots). For a comparison, we also show the nine galaxies with both Lya and Ha emission in 
ISong et ^ j20l4) with blue filled dots (UV continuum slopes are not provided for individual galaxies in ISong et alj2014l) . The Lya escape fraction has been 
calculated by using the SED-derived dust attenuation and assuming that Es{B — V) = Eg(B — V). This assumption is the main factor affecting the uncertainties 
of the results, not only in this work, but also in previous results in the literature since the relation between Es{B — V) and Eg{B — V) at z ~ 2 has not been 
accurately established yet. The straight lines are linear fits to our points. We also include upper limits for the HAEs undetected in Lya. It can be seen that the 
Lya escape fraction decreases with increasing dust attenuation, redder UV continuum slopes and higher stellar masses and SFRs. 


when galaxies are averaged over different bins of Ha-derived SFR 
or stellar mass. This result shows the robustness of the Ha emis¬ 
sion as tracer of SFR at z ~ 2, which adds to the unbiased and 
well-understood selection functi on of the Ha NB method to find 
SF galaxies at different redshifts dSobral et alJl2012l) . 

Previous works have also analysed the accu racy of Ha emis¬ 
sio n to recover SFR at high redshift, for example [Erb et alj ( l2006h 
and iReddv et an( l2010h . They correct from dust attenuation derived 
from SED fitting, Es{B~V), and assuming Es{B~V)= Eg{B — V), 
where — V) is the reddening for nebular emission. They actu¬ 
ally obtained tha t the traditionally em ployed relation Es{B — V) = 
0.4 X Eg{B — V) dCalzetti et alj|200^ produces Ha-SFR that over¬ 
predict those de rived from the X-ra y and dust-corrected UV emis¬ 
sions. Recently. ISteidel et alj d2014h apply a relation between Ana 
and Es{B — V) that depends on the value of Es{B — V). The main 
advantage of our dust correction method for Ha emission is that 
we do not require any assumptions about the relationship between 
Es{B — V) and Eg{B — V), which is still until debate and have been 
suggested to be redshift-dependent dKashino et al.ll201^ . Instead, 
we use a relation between dust attenuation and stellar mass that 
has been reported to be valid at least up to z ~ 1.5 for HAEs 
JSobral et alj|2012l) . 

We also compare in the right panel of Figure|^the Ha-derived 
total SER with the total SER obtained with direct Herschel detec¬ 
tions for the 9 PACS/SPIRE-detected HAEs. As it can be seen, the 
Ha emission, after the dust correction has been included, is not able 
to recover the more accurate SFR derived with PACS detections. 


This is because, at z ~ 2, Herschel only detects the most extreme 
sources for which our average dust correction factor applied to the 
Ha luminosity are not high enough (due to very high internal ob¬ 
scuration) to reproduce the more accurate Herschel-deiived SFR. 
This also happens to the dust correction factors derived from the 
UV continuum s lopes in high-red s hift Herschel-selected galaxies 
(see for example loteo et al.l 2013al : lR'odighiero et al.ll2014l) . Stack¬ 
ing as a function of the Ha-derived total SFR we only recover one 
> 3 ct stacked detection, also represented in the right panel of Figure 
|9] The total SFR derived from the PACS stacked flux is higher (in 
about ~ 0.3 dex) with respect to the Ha-determination. This might 
indicate that the dust-correction factor used to correct the Ha emis¬ 
sion is more uncertain for the most massive HAEs with the highest 
SERs, although deeper EIR data would be needed to confirm this 
trend to more normal SE galaxies with lower SER. 


6 SFR VERSUS STELLAR MASS RELATION AND ITS 
UNCERTAINTIES AT Z ~ 2 

6.1 The main sequence for HAEs at z ~ 2 

Most previous works agree that there is a relation between SER 
and stellar mass, commonly referred to as the main sequence ( MS), 
where normal SE galaxies are located (see lSpeagle et al .12014 for a 
recent compilation). The MS has been reported to exist from the lo¬ 
cal universe u p to high-redshift an d to be relatively independent of 
environment ( IKovama et al.ll2(H3h . However, less is known about 
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Figure 9. Compaiison between the SFR derived from Ha and rest-frame UV luminosities for our studied HAEs at z ~ 2. Left panel: results without dust 
correction. Right panel: results w hen dust correction is included. The Ha luminosities are correcte d by dust attenuatio n by using the local relation between 
dust attenuation and stellar' mass jOam & Beslll201ol) . s uggested to be valid at least up to z ~ 1.5 iSobral et al.ll201 Jl . The rest-frame UV luminosities are 
corrected by using the dust attenuation derived with the iHeinis et alj 1201 3h relation. In both panels, the filled red squares represent the median SFRuv in 
different bins of SFRna and filled cyan squares represent the median SFRuv in different bins of stellar mass. The solid lines are the one-to-one relations and 
the dashed lines are deviations of ±0.3 dex. PACS-detected and stacked HAEs are represented in the right panel by filled green and orange dots, respectively. 
The enor bars on the right panel shows a lower limit o n the uncerta i nties of the results. Eor the dust correction applied to Ha luminosity, this uncertainty 
comes from the scatter in the best-fitted relation found in ISobral et al.l I 20 T 2 I) . It should be pointed out that the uncertainties affecting to the dust con'ections in 
Ha are much lower than those affecting to the rest-frame UV continuum, since the latter is much more affected by dust. 


the values of its zero point and slope at a given redshift. One of the 
main reasons is the lack of FIR detections for a representative pop¬ 
ulation of SF galaxies at each redshift, that prevents accurate deter¬ 
minations of the total SFR. Ev en with the deepe st FIR surveys car¬ 
ried out so far (see for example felbaz et alj20I ih . only a very small 
fraction o f the galaxies are detected in the FIR, mostly at the highest 
redshifts iOteo et alj|2013al[2014t) . Furthermore, it is claimed that 
most high-redshift FIR-detected galaxies are not normal SF galax¬ 
ies, but likely have a starburst nature and are preferentially located 
above the MS teodighiero et alj|201 ll : lUee et alj|^013l : lOteo et alj 
l2013al . l2014h . Therefore, even with Herschel-sAecX&A galaxies it is 
not possible to determine the slope and zero point of the MS at high 
redshift. 

Now that we have evidence from ^that the Ha emission is 
a good tracer of star formation at z 2.23, that Ha provides a 
clean, well-understood, SFR-selected sample (not biased to either 
just blue or just red galaxies), that HAEs are an excellent repre¬ 
sentation of the whole population of SF galaxies at z ~ 2(0, and 
that their stellar masses cover a wider range than many previous 
studies at z ~ 2, we can attempt to study the relation between SFR 
and stellar mass and its uncertainties by using our sample of HAEs. 
This is shown as the orange points in Figure [TO] with the best-fit 
MS indicated by the red line. 

The main caveat of the analysis of the MS with our sample 
of HAEs is that they are s elected down to a fi xed SERHa-imcon- 
(~3M0yr^* and EW cut. ISobral et al.lUoi3h . Furthermore, the 
SFRea-corr of HAEs is obtained by using a dust correction factor 
which depends on the relation between dust attenuation and stellar 


mass. All of this means that, for a given stellar mass, only HAEs 
with SERna-corr above a certain limit (that is dependent on stel¬ 
lar mass) are included in the sample. This might produce a flat¬ 
tening of the MS with respect to previous works based on other 
selection functions. In order to explore this in more detail, we have 
obtained the limiting SFRea-corr in different bins of stellar mass. 
This is represented in Figure [TO] with the open black squares and 
the dashed curve. As expected, the limiting curve tracks along the 
locus of the bottom of the points corresponding to HAEs. As an 
example, our Ha selection would no t be able to detect galaxies on 
the MS relation of lOaddi et alj ( l200^ solid black line on FigurefTOt 
that have stellar masses log [M^/Mq) ~ 9.5 (but note t hat dust cor¬ 
rection factors used to recover total SFR are different in lPaddi et alj 
l l2007l) and this work). If our Ha observations would have selected 
galaxies down to a lower SFRea-uncorri we could have been able 
to select galaxies with lower SFRna-coiT at that stellar mass and, 
consequently, we might have obtained a steeper MS. This reflects 
the effect of our Ha selection on the slope of the MS at z ~ 2 and 
should be taken into account in any comparison with the results 
presented in the literature for galaxies selected in a different way. 
However, as it will be shown in the next Section, when we use the 
same methods that in other works to determine the SFRna-coiT and 
stellar mass, the MSs are in very good agreement. 

It is clear from Figure [T^ that an MS for HAEs exists at 
logiMif/Mo) > 9.25: more massive HAEs have higher SERs on 
average. However, the scatter is significant and it increases with 
stellar mass. It can be also seen that the MS is flatter for more 
massive galaxies. Actually the relation becomes almost flat at 
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Figure 10. Relation between the SFR and stellar mass for our HAEs at 
z ~ 2.23 (orange dots). Light blue squares represent the median SFR in dif¬ 
ferent bins of stellar mass. The SFR of HAEs is derived with the lKennicud 
h998l) calibration and the Hct luminosity con'ected from dust attenuation 
with the local dust-mass relation. The location of Herschel-^oiQciQd and 
stacked HAEs is indicated with the orange diamonds and the brown filled 
dot, respectively, and their total SFR are SFRtotal — SFRuv + SFRir. For 
compaiison, we also show a sample of stacked sBzK galaxies with filled 
grey dots. LAEs are represented with green dots and their SF R are derived 
from the dust-corrected rest-frame UV luminosity with the iHeinis et alj 
J2013t) relation. Previously reported MS at z ~ 2 are over-plotted, as indi¬ 
cated in the bottom-right legend. For all the galaxies, the stellar masses have 
been calculated from SED fitting with BC03 templates including emission 
lines. The open black squares and dashed curve represent the limiting dust- 
conected SFR as a function of mass for our Ha survey (see text for more 
details). It can be seen from this figure that the HAEs follow an MS with a 
lower slope than the MS at z ~ 2 reported in previous works. Furthermore, 
our results indicate that there is a flattening of the MS at the highest stel¬ 
lar masses, as previously reported in previous works from different analysis 
based on different samples of SF galaxies. 


log(M*/M0) > 10.2, where the width of the MS is also higher. 
Since HAEs are selected by star formation, they are truly SF galax¬ 
ies: all HAEs with log(M*/M0) ^ 10.0 would also have been se¬ 
lected as SF galaxies according to the sBzK criterion. Therefore, 
the flattening is not caused by the presence of massive quiescent 
galaxies without star formation. Thi s flattening of the M S towards 
massive galaxies can be also seen in iHeinis et alJ l l2014l) at z 1.5 
in their work about stacking analysis in Herschel bands, although 
the flattening found in this work is much more significant (see also 
IOteoll2014llWhitaker et al.ll20l4 or lSchreiber et al.l[2014h . By fit¬ 
ting a linear relation for galaxies with log(M*/Af0) ^ 10.25 in 
the form log (SERna-corr) = a-\-b x log(M,/M0) we find a = 
—3.65 ±0.15 and b = 0.52 ±0.02. This MS is in very good agree¬ 
ment with the results shown in IZahid et ^ (1201 2h in a sample of 
spectroscopically confirmed Ha emitters over a similar stellar mass 
range than our sample. This is the most comparable sample to ours 
that can be found in t he literature in terms of both Ha emission 
and stellar mass range. [Kashino et ^ ( 1201 3l) obtained a higher MS 
slope with a spectroscopically confirmed sample of galaxies with 


Ha emission, but their galaxies are at a slightly lower median red- 
shift than HAEs in this work and their sample is restricted to more 
massive galaxies, mainly with \og{M^/MQ) >10. Our MS has 
a lower slope than the classical MS found in IPaddi et"^ ( l2007h 
(D07), probably produced by our selection effect, although differ¬ 
ent dust-correction factors and recipes for stellar mass determina¬ 
tion might play a significant role. This will be discussed in more 
detail in 36.21 

PACS-detected HAEs are all well above the HAE MS, as hap¬ 
pens for PACS-de tected LBGs and sBzK galaxies dOteo et alj2014l : 
iRodighiero et aljEzOlA . This indicates that they have an SB nature 
and do not belong to the population of normal SF galaxies at z ~ 2. 
Actually, our PACS-detected sB zK galaxies have even higher SFRs 
than those in[Oteo^^h ( 201^ due to the shallower PACS data 
in COSMOS dLutz et alJl2oTil) . When stacking as a function of 
stellar mass, we only recover one detection for HAEs, also rep¬ 
resented in Figure [Tol The stacked point is in agreement with the 
D07 MS and it is located slightly above the MS for HAEs found 
in this work. This is compatible with the previous result that the 
SFR derived from stacking is slightly higher than the derived from 
Ha in the most massive HAEs. For comparison, we also show the 
stacked points corresponding to the sample of sBzK galaxies. Only 
the most massive sBzK galaxies are detected through stacking. The 
stacked points for sBzK galaxies agree very well with the extrapola¬ 
tion of the MS for log {Mi,/MQ) ^ 10.25 HAEs towards higher stel¬ 
lar masses. However, the most massive stacked sBzK galaxies have 
higher SFR than the most massive HAEs. This difference might be 
due to the fact that the dust correction factor employed to recover 
the total SFR might not be accurate for the most massive galaxies. 
However, it could also be as a consequence of a different nature be¬ 
tween HAEs and sBzK galaxies for the most massive objects. The 
latter would be supported by the shape of their UV-to-NIR SEDs 
(see Figure |7][, since massive HAEs are bluer than massive sBzK 
galaxies. 

As already obtained in ^ galaxies selected through the Lya 
technique tend to be less massive than HAEs. Furthermore, the 
Lya selection allows to probe down to lower dust-corrected SFRs. 
This can be used to extend the MS obtained for HAEs down to 
\og{Mt/M q) ~ 9. This is represented by the blue filled squared 
in Figure [Tol As in ^ the dust-corrected SF R for LAEs has bee n 
obtained by assuming the IRX-p relation of iHeinis et alj ( l2013h . 
Interestingly, and despite the effect of our Ha selection on the MS 
slope, the point for LAEs agrees very well with the extrapolation 
of the HAE-MS to lower stellar masses. 


6.2 Uncertainties in the SFR-mass relation 

In order to study the uncertainties involving the determination of 
the slope and zero-point of the HAE MS at z ~ 2 and try to explain 
in more detail the differences found with previous works, we com¬ 
pare now the SFR-mass relation when considering different recipes 
used in the literature to obtain stellar mass and total SFR. We fo¬ 
cus first on the dependance of the HAE MS upon the recipe used 
for determining stellar mass. This is shown in Figure [TT] where all 
the SFRs have been derived by c orrecting the rest-fram e UV lumi¬ 
nosity with the IRX-P relation of iMeurer et alj ( Il999h for the sake 
of homogeneity with previous works that will be referred to in this 
Section. Therefore, the main difference would be only the determi¬ 
nation of stellar mass. We include in Figure [TT] the MS when the 
stellar mass of HAEs is obtained through our method of SED fits 
with BC03 templates associated to exponentially declining SFHs 
and the inclusion of emission lines (blue open squares, see 
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Figure 11. Relation between dust-corrected total SFR and stellai* mass for 
our HAEs when different recipes to derive stellar mass are used: SED fits 
with BCQ3 temp l ates a nd inclusion of emission lines (blue open squai'es) 
and lOaddi et~^ i2QQ4) calibration between stellai' mass and rest-frame 
optical color (orange filled dots). In all cases we derive dust-corrected 
total S FR from the r est-fra me UV luminosity corrected with the rela¬ 
tion of iMeurer et alj jl999l) IRX-p. For compai'ison, we also show the 
points associated to stacked sBzK galaxies (grey open dots) and HAEs (red 
open dot) in PAC S-160 /.im whose stellar masses are obtained from the 
lOaddi et ai] j2004h calibration. Inset panels show the relation between stel¬ 
lar masses when different approaches are consider ed for their calcula tion: 
SED-derived mass against masses obtained with the lPaddi et al.ll2Q04) cal- 
ibration (top-left panel) and SED-derived median versus best-fitted masses 
(bottom-right panel). In both inset panels, black line is the one-to-one re¬ 
lation, and the red line is the linear fit to the points. This figure indicates 
that the definition of the MS (slope and normalization) depends on the way 
stellar mass is calculated. 


and also (orange dots) that obtained using the calibration between 
BzK phot ometry and stellar mas s reported in iDaddi et al.l j2QQ4t) 
(see also iRodighiero et al.l|2Q14l) . It can be clearly seen that the 
slope of the HAE MS depends upon the method used for deriv¬ 
ing stellar mass, with the slope being lower when stellar masses are 
derived with SED fits with inclusion of emission lines. The top-left 
inset plot of Figure fTTI compares the stellar masses derived with 
BC03 templat es including the effe ct of emission lines with the de¬ 
rived with the [Paddi et"^ ^2004) calibration. It can be seen that 
they are correlated but do not follow the one-to-one relation, ex¬ 
plaining the difference in the MS. 

When using the same dust correction factors and recipe for 
stellar mass determination than in D07, we find a MS slope closer 
to their result (cf. orange and black lines in Figure im . The still 
slightly lower MS slope might be explained by the Ha selection 
effect (see We also plot in Figure [TT] the stacked points for 

HAEs and sBzK g alaxies when stellar masses are obtained as in 
IPaddi et'H] ll2004l). It can be seen that they all follow quite well 


Paddi et alj j2007h MS. indicatin g that if stellar mas ses are cal¬ 
led in the same manner then thelPaddi et alj (l2Q07h MS is re- 


the 

culated _ 

cover ed from different SFR estimators (see also lRodighiero et alj 
l2014h . We have also compared the best-fitted stellar masses with 


the median values as derived by LePhare (bottom-right panel in 
Figure [YTT l. However, no significant differences are seen between 
both results that might significantly alter the slope of the HAE MS 
at z ~ 2. 

Eigure [T^ shows the impact of different dust correction fac¬ 
tors in the definition of the HAE MS at z ~ 2. With the aim of 
avoiding any uncertainty coming from the determination of stel¬ 
lar mass we represent the SFR-mass relation by using the rest- 
frame K band luminosi ty. This luminosity is a proxy of stellar 
mass llProrv et al.l 12004)) and is a direct observable quantity that 
does not need any assumption or calibration in its determination. 
Three different methods for deriving the total SER have been em¬ 
ployed: 1) dust-corrected Ha luminosity with the local dust-mass 
relation (red points); 2) d ust corrected rest-f rame UV luminosity 
with the IRX-p relation of lHeinis et alJ ( l2013l) (orange filled dots); 
3) dust corrected res t -frame UV luminosity with the IRX-p rela¬ 
tion of iMeurer et al.l d 19991) (blue filled dots). Note that the only 
difference in Figure [T^ are the methods adopted for deriving the 
total SFR. It can be seen that the MS is strongly affected by the dif¬ 
ferent dust corrections, with the slope being steeper with the dust 
correction based on the IMeurer et al.l i ll 9991) relation. The MS de¬ 
rived from the Ha emission assuming the local dust-mass relation 
and from the rest-frame UV w ith the dust correction based on the 
relation of lHeinis et alj (1201 3l) are similar due to the agreement be¬ 
tween both SFR indicators (see Figure|9]l. For a reference, we show 
in the inset panel of FigurefT^the relation between stellar mass and 
rest-frame K band luminosity for the two different assumptions for 
stellar mass determination used in the discussion above. 

Summarizing, our results show that, despite there is evidence 
that the MS exists, its slope and zero points at each redshift are 
challenging to determine. Their values depend on the recipes em¬ 
ployed to obtain stellar mass and total SFR. For given assumptions 
of these, consistent values can be determined provided that repre¬ 
sentative samples of sources such as NB-selected HAEs or sBzK 
galaxies are used for their characterization. 


7 CONCLUSIONS 

In this work, we have carried out a multi-wavelength analysis (from 
the rest-frame UV to the FIR) of the spectral energy distribution 
(SED) of narrow-band (NB) selected, star-forming (SF) Ha emit¬ 
ters (HAEs) with the aim of analysing their physical properties and 
their importance for galaxy formation and evolution. We have also 
compared their physical properties with those derived for other 
classical populations of SE galaxies at their same redshift: NB- 
selected Lya emitters (LAEs), Lyman-break galaxies (LBGs), and 
sBzK galaxies. The main conclusions of our work are the following: 

(i) The HAE selection recovers the full diversity of SE galaxies 
at z ~ 2. Coupled with the simple and well-understood HAE selec¬ 
tion (selecting SF galaxies down to a given dust-uncorrected SFR) 
which can be self-consistently applied at multiple redshift slices 
from the local Universe up to z ~ 2.5, HAEs represent an excellent 
sample to study galaxy evolution. 

(ii) At the depth of the COSMOS data, only about 30% of sBzK 
galaxies can be selected with the drop-out technique, whereas 95% 
of LBGs can be selected with the sBzK criterion. Only 4.5% of 
HAEs can be selected as LAEs. These results indicate that the 
Lyman-break and Lya selections miss a relevant percentage of the 
SF population at the peak of cosmic star formation. LBGs and 
LAEs are biased towards blue, less massive galaxies. Although the 
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Figure 12. Relation between total SFR and the rest-frame K band lumi¬ 
nosity for our sample of HAEs at z ~ 2.23. Different estimators of the total 
SFR have been used: dust-corrected Ha luminosity with the local dust-mass 
relation (red dots) and dust-co n'ected UV rest-fram e luminosity with the 
iHeinis et al ] j2Q13h (blue) and iMeurer et al.l jl999h IRX-p (orange) rela¬ 
tions. The inset panel shows the relation between the stellar mass and the 
rest-frame K band luminosity for two different recipes for stellai' mass: SED 
fits with BC03 tem plates (green) and calibration with z and K magnitudes 
jPaddi et a 1 J 2004 (blue). The black line would correspond to mass-to-light 
ratio equal to unity. Green and blue straight lines are fitted to the green and 
blue points, respectively. This figure indicates that dust-correction factors 
have a strong influence on the definition of the MS at z ~ 2, both in the 
slope and zero point. 


precise numerical results will depend on the depth of the observa¬ 
tions considered in the analysis, these results are important to in¬ 
terpret results at higher redshifts where only Lya or Lyman-break 
selections can be applied. On the other hand, most LBGs and LAEs 
can be selected through the sBzK criterion. Only the bluest LAEs 
and LBGs would be missed, although this sample represents a very 
low percentage of the SE population at z ~ 2 and their exclusion 
would not affect significantly conclusions for galaxy evolution. 

(iii) There is a significant percentage of LAEs that are not de¬ 
tected in optical and near-IR broad band filters even with the deep 
photometry used in this work. The non detection indicates that 
these galaxies have a faint continuum but strong emission lines. 
Due to their non detection, SED fits cannot be carried out for these 
galaxies and, consequently, their properties and, most importantly, 
their contribution to galaxy evolution studies, are unknown. 

(iv) Although the Lya criterion preferentially selects SE galax¬ 
ies with low dust attenuation and low stellar mass (likely due to the 
resonant nature of the Lya emission), there is also a small percent¬ 
age of red and massive LAEs in our sample, even redder than any 
LBG in our LBG sample. This indicates that Lya is also able to es¬ 
cape from dusty and massive galaxies, in agreement with previous 
work reporting Lya emission in sub-mm and EIR selected samples. 
Since the number of red LAEs is low compared with the total pop¬ 
ulation, surveys over wider areas are needed to study in detail the 


properties of these galaxies and shed light on how Lya can escape 
from these systems. 

(v) The median SEDs of the galaxies studied in this work reveal 
that the Lya emission is strong enough to affect the broad-band U 
photometry only in LAEs. This result indicates that Lya emission 
is not strong in LBGs, HAEs, or sBzK galaxies on average, high¬ 
lighting the low probability of finding Lya emission in a general 
sample of selected SE galaxies. 

(vi) Only 4.5% of HAEs show detectable Lya emission imply¬ 
ing low Lya escape fraction at z ~ 2 in agreement with previous 
results. Additionally, we find that the Lya escape fraction (/esc) de¬ 
creases with increasing SED-derived dust attenuation, the UV con¬ 
tinuum slope, stellar mass and SER. This suggests that Lya pref¬ 
erentially escapes from blue galaxies with low dust attenuation, al¬ 
though a population of red LAEs is also present indicating that dust 
and Lya are not mutually exclusive. 

(vii) By using completely different and independent methods to 
recover the total SER, we obtain that the Ha emission is an excel¬ 
lent star formation tracer at z 2 with deviations typically lower 
than 0.3 dex for individual galaxies. These deviations are close to 
zero when averaging the sample in stellar mass or SER bins. 

(viii) By using the Ha-derived SFR we study the relation be¬ 
tween SFR and stellar mass for our HAEs. We find a main-sequence 
(MS) of star form ation, but with a slope lower than the classical 
IPaddi et"^ l l2007h relation. By fitting a linear function in the form: 
log SFR = a + bx log (M/Mq), we obtain: a = —3.65 ± 0.15 and 
b = 0.52 ±0.02. We show that, in part, the lower slope with respect 
to previous works might be due to the different selection criteria. 
However, exploring the uncertainties in the slope and zero point of 
the HAE MS, we find that they are very sensitive to both the dust 
correction factors adopted to recover the total SFR and the way the 
stellar masses are determined. This largely explains the difference 
with previous works and represent the main uncertainty in the def¬ 
inition of MS at high redshift. This might apply to any sample of 
SE galaxies. Using t he same recipes for s tellar mass calculation and 
dust correction as in lPaddi et alj (|2007|) we find consistent results 
for the MS at z ~ 2 for our HAEs and also for sBzK galaxies whose 
total SFR are obtained with a stacking analysis in Herschel bands. 
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